Airborne measurements of C0 2 column absorption and 
range using a pulsed direct-detection integrated 
path differential absorption lidar 


James B. Abshire, 1 * Haris Riris , 1 Clark J. Weaver , 2 Jianping Mao , 2 Graham R. Allan , 3 
William E. Hasselbrack , 3 and Edward V. Browell 4 

^ASA Goddard Space Flight Center, Greenbelt, Maryland 20771, USA 
2 Earth System Science Interdisciplinary Center, University of Maryland, College Park, Maryland 20740, USA 
3 Sigma Space Corporation, Lanham, Maryland 20706, USA 
4 STARSS-II Affiliate, NASA Langley Research Center, Hampton, Virginia 23681, USA 
Corresponding author: James.B.Abshire@nasa.gov 

Received 31 January 2013; revised 25 March 2013; accepted 11 April 2013; 
posted 6 May 2013 (Doc. ID 184647); published 21 June 2013 


We report on airborne C0 2 column absorption measurements made in 2009 with a pulsed direct- 
detection lidar operating at 1572.33 nm and utilizing the integrated path differential absorption technique. 
We demonstrated these at different altitudes from an aircraft in July and August in flights over four 
locations in the central and eastern United States. The results show clear C0 2 line shape and absorption 
signals, which follow the expected changes with aircraft altitude from 3 to 13 km. The lidar measurement 
statistics were also calculated for each flight as a function of altitude. The optical depth varied nearly 
linearly with altitude, consistent with calculations based on atmospheric models. The scatter in the optical 
depth measurements varied with aircraft altitude as expected, and the median measurement precisions 
for the column varied from 0.9 to 1.2 ppm. The altitude range with the lowest scatter was 8-10 km, and 
the majority of measurements for the column within it had precisions between 0.2 and 0.9 ppm. 

OCIS codes : (280.3640) Lidar; (300.6360) Spectroscopy, laser; (010.3640) Lidar; (140.3510) Lasers, 

fiber. 
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1. Introduction 

Although increasing atmospheric C0 2 is widely ac- 
cepted as the largest anthropogenic factor causing 
climate change, there is considerable uncertainty 

about its global budget. To be useful in reducing 
uncertainties about carbon sources and sinks, 

space-based atmospheric C0 2 measurements need 
-0.3% precision on spatial scales of <10 x 10 deg 
[1-3]. Several groups have analyzed space missions 
using passive spectrometers [4-6], and the GOSAT 
mission [7] is now making global C0 2 measurements 
from space using a Fourier-transform infrared spec- 
trometer and reflected sunlight. However, sun angles 
restrict its measurements to the daytime and 


midlatitudes. An inherent error source with space- 
based passive spectrometers is variability in path 
length caused by changes in surface elevation and 
by optical scattering from aerosols and thin clouds in 
the illumination or observation paths [8,9]. Scatter- 
ing in either of these paths modifies the optical path 
length and thus the total C0 2 absorption viewed by 
these spectrometers, and so it can cause large 
retrieval errors even for thin cirrus clouds [10]. 

To overcome these limitations, the US National 
Research Council’s 2007 Decadal Survey recom- 
mended a new space-based C0 2 measuring mission 
called ASCENDS [11]. Its goals are to produce global 
atmospheric C0 2 measurements using the laser ab- 
sorption spectroscopy approach. The mission’s goals 
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are to quantify the global spatial distribution of 
atmospheric C0 2 with 1-2 ppm accuracy and to 
quantify the global spatial distribution of terrestrial 
and oceanic sources and sinks of C0 2 with approxi- 
mately monthly time resolution. The ASCENDS 
approach offers continuous measurements over the 
cloud-free oceans, at low sun angles and in darkness, 
which are major advantages over passive sensors. 
The ASCENDS mission organizers held a workshop 
in 2008 to better define the science and measurement 
needs and plans for future work [ 12 ] . The European 
Space Agency (ESA) has also conducted mission 
definition studies for a similar space mission called 
A-SCOPE [3,13], and their lidar sensitivity and 
spectroscopic analyses have recently been publi- 
shed [14,15] . 

2. C0 2 Line Choice and Lidar Approach 

The integrated path differential absorption (IPDA) 
lidar technique is widely used for open-path laser 
absorption spectroscopy measurements [16,17] . 
Typically two laser wavelengths are used, and the 
target is illuminated with the laser alternatively 
tuned onto and off the gas absorption line. The 
path-integrated gas absorption attenuates the on- 
line laser energy relative to the off-line wavelength. 
By measuring the optical depth of the gas absorption 
line, along with path length, and by knowing the 
difference in gas absorption cross sections, one can 
solve for the path-integrated gas number density. 
Several groups have measured atmospheric C0 2 
using IPDA lidar, and summaries of prior work have 
recently been published [ 18 , 19 ]. 

The motivation for developing our airborne lidar 
[ 18 , 20- 23] is to demonstrate a dual-wavelength 
pulsed IPDA approach as a candidate for the 
ASCENDS mission. The approach uses two tunable 
pulsed laser transmitters allowing simultaneous 
measurement of the absorption from a C0 2 absorp- 
tion line in the 1570 nm band, 0 2 extinction in the 
oxygen A-band, and surface height and atmospheric 
backscatter in the same path. A tunable laser is 
stepped in wavelength across a single C0 2 line for 
the C0 2 column measurement, while simultaneously 
a laser is stepped across a line doublet near 765 nm in 
the oxygen A-band for an atmospheric pressure meas- 
urement [24,25] . Both lasers are pulsed at a ~8 kHz 
rate, and the two absorption line regions are sampled 
in wavelength steps at typically ~1 kHz rate. The la- 
ser transmitters utilize tunable diode lasers followed 
by laser amplifiers. The direct-detection receiver mea- 
sures the time-resolved laser backscatter from the 
atmosphere along with the energies of the laser 
echoes from the surface. After suitable averaging, 
the gas extinction and column densities for the C0 2 
and 0 2 gases are estimated from the sampled wave- 
lengths of the surface-reflected line shapes. 

We use a single line in the 1570 nm band for the C0 2 
measurement [9]. This vibration-rotation band of 
C0 2 has an appropriate range of absorption that 
provides good sensitivity to the surface echo signal 


and to variation in C0 2 in the lower troposphere. It also 
has minimal interference from other atmospheric 
species such as H 2 0, and has several temperature- 
insensitive lines. The shorter wavelength lines in 
the R branch are a better match to available laser 
and detector technologies. The centermost line of 
the R branch at 1572.335 nm has been analyzed and 
is attractive for C0 2 measurements [26]. It has a 
low temperature sensitivity, particularly to tempera- 
ture changes in the lower atmosphere. 

Our approach samples the C0 2 line shape at multi- 
ple wavelengths. Nominally, eight are planned for 
space, while we have used both 20 and 30 for 
airborne measurements. Sampling at multiple wave- 
lengths across the absorption line provides several 
capabilities and supports analysis using at least 
two vertical atmospheric weighting functions, which 
peak separately at lower and upper altitudes [26]. 
This approach also provides significant additional in- 
formation that allows for solving for instrument 
parameters, such as baseline tilts, wavelength off- 
sets, and wavelength-dependent lidar transmissions, 
which allows optimization of the measurements in 
subsequent analysis. To date we have used this for 
solving for wavelength offsets via a line-fitting proc- 
ess. We have also used it to characterize the lidar’s 
response as a function of wavelength, which allows 
solving for and reducing the impacts of any residual 
wavelength variability. Oversampling the line in air- 
borne lidar compared to space also allows postmea- 
surement analysis to vary the number of samples 
used and to help determine the optimum number. 

Using pulsed lasers and receiver processing to 
time (and range) resolve the laser backscatter pro- 
files allows postdetection signal processing to isolate 
the laser echo signals from the surface, and to reject 
backscatter from the atmosphere that arrives earlier. 
Hence it allows isolating the full column measure- 
ment from bias errors caused by atmospheric scatter- 
ing [8,9] . Isolating the surface echo with a time gate 
in receiver processing also substantially improves 
the receiver’s signal-to-noise ratio (SNR) by reducing 
the amount of noise from the detector and the solar 
background. For space missions, the photon effi- 
ciency of the lidar measurement technique is quite 
important. The photon efficiencies for three IPDA 
lidar approaches have recently been analyzed and 
compared [27] and the results showed that the pulse 
technique was quite efficient. 

3. Airborne Lidar Description 

In 2008 we adapted our ground-based lidar [20] for 
airborne use on the NASA Glenn Lear-25 aircraft, 
[28] shown in Fig. 1. The sensor head was mounted 
above the aircraft’s nadir window. In initial flights 
[ 18 ] , the etalon fringes from the aircraft’s nadir cam- 
era window significantly limited the measurement 
precision. For the 2009 flights a custom pair of 
wedged and antireflection coated optical windows 
were used as a replacement, and they eliminated this 
issue. 
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Fig. 1. (Left) NASA Glenn Lear-25 aircraft. The nadir window assembly is just below the NASA logo. Photographs of the lidar installed on 
the aircraft showing the sensor head assembly (middle) and the dual aircraft racks (right). 


The airborne lidar transmitter operates in pulsed 
mode by using an acousto-optic modulator (AOM) be- 
tween the diode laser and the fiber amplifier. The 
receiver measures the time-resolved sum of 450 
photon-counting lidar backscatter profiles each sec- 
ond, for each of the 20 laser wavelengths, with suffi- 
cient recording duration to include the surface echo 
signal. The airborne lidar specifications are listed in 
Table 1. 

The approach for sampling the C0 2 line and lidar 
block diagram is shown in Fig. 2. The lidar signal 
source is a distributed feedback (DFB) laser diode, 
which is operated near 1572.33 nm by controlling 
its temperature and current. A ramp from a signal 
generator was used to sweep the current to the diode 
laser, and hence its output wavelength. The diode’s 
CW output is then gated into pulses by an AOM to 
an erbium-doped fiber amplifier (EDFA). For these 
measurements, 20 wavelength samples are used 


Table 1. 2009 Pulsed Airborne C0 2 Lidar Parameters 


C0 2 line center wavelength 

1572.33 nm (typically) 

Laser min & max wavelengths 

1572.29 nm, 1572.39 nm 

Laser wavelength steps 

20 (these flights) 

across line 

Laser wavelength change/step 

-5 pm 

Laser peak power, pulse width, 

25 W, 1 ps, 25 pJ 

energy 

Laser divergence angle 

100 prad (these flights) 

Seed laser diode type 

DFB: Fitel FOL15DCWD 

Laser pulse modulator (AOM) 

NEOS Model: 26035-2-155 

Fiber coupled C0 2 cell 

80 cm path, -200 Torr 

Fiber laser amplifier (EDFA) 

pressure 

IPG EAR-10K-1571-LP-SF 

Laser line scan rate 

450 Hz 

Laser linewidth per step 

-15 MHz 

Receiver telescope 

Cassegrain, f / 10 (Vixen) 

Telescope diameter 

20 cm 

Receiver FOV diameter 

200 prad 

Receiver optical bandwidth 

800 pm FWHM 

Receiver optics transmission 

0.64 

Detector PMT type 

Hamamatsu H10330A-75 

Detector quantum efficiency 

2% (this device) 

Detector dark count rate 

-500 kHz 

Receiver signal processing 

Photon counting/ 

Histogram time bin width 

histogramming 
8 ns 

Receiver integration time 

0.9 s per readout 

Recording duty cycle 

45% (0.9 s every 2 s) 


across the C0 2 line. A small percentage of the CW 
seed laser output is split off and directed through 
a fiber-coupled C0 2 absorption cell and to a detector. 
The C0 2 cell serves as a monitor for the center wave- 
length of the sweep. The nominal static transmitted 
laser pulse wavelengths are measured in a calibra- 
tion procedure using a commercial wavemeter. Sub- 
sequent testing showed some curvature in the 
scanning dynamic ramp signal, so we used a more 
accurate quadratic functional model for the laser 
wavelength versus pulse sequence number in the 
data analysis. 

The laser output is a sequence of 1 ps wide laser 
pulses that occur at a 10 kHz rate. The collimated 
transmitted laser signal exits toward nadir though 
the window. The laser backscatter is collected by 
the receiver’s 20 cm diameter Cassegrain telescope, 
which views nadir through another window in a bi- 
static configuration. A multimode fiber is used to cou- 
ple the optical signal from the telescope focal plane to 
the receiver optics. After passing through an optical 
bandpass filter, the signal is focused onto a photo- 
multiplier tube (PMT) detector. The PMT used in 
the 2009 flights had a single photon-counting effi- 
ciency of ~2% near 1570 nm. The electrical pulse out- 
put from the PMT was amplified and passed through 
a threshold detector. 

The pulses from the receiver’s PMT detector and 
discriminator are binned and accumulated by the 
multichannel scalar (MCS). The start time of the 
MCS sweep is synchronized with the trigger for 
the first laser pulse of the 20 and hence the start 
of the pulsed wavelength sweep. Each MCS sweep 
contains a time-resolved histogram of PMT pulse 
counts versus time (i.e., the laser backscatter pro- 
files) for all 20 pulses. At the end of 0.9 s, each 
MCS bin contains the total receiver counts for the 
450 laser sweeps. The receiver histogram record is 
then read and stored. The laser trigger and data ac- 
quisition are synchronized to timing markers from 
the GPS receiver, and data was stored every other 
second. The computer also digitizes other signals, in- 
cluding those from thermocouples distributed across 
the sensor head and electronic rack, the inertial guid- 
ance system output from the aircraft, and the GPS 
position and time. A nadir viewing video camera also 
records the visible image though the nadir window. 
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Fig. 2. (Left) C0 2 line sampling approach for airborne lidar. (Right) C0 2 lidar block diagram. 


4. 2009 Airborne Campaigns 

During mid-July though mid-August 2009, we flew 
an airborne campaign with a series of flights, each 
about 2.5 h long. We measured the shapes of the 
atmospheric C0 2 absorption line at stepped altitudes 
from 3 to 13 km over a variety of surfaces in 
Nebraska, the Department of Energy’s ARM site 
in Lamont, Oklahoma, central Illinois, and near 
and over the Chesapeake Bay in North Carolina 
and Virginia’s Eastern Shore Peninsula. These 
flights allowed testing and recording performance 
under different measurement conditions. 

Strong laser signals and clear ground echo pulses 
and C0 2 line shapes were observed at all altitudes on 
most flights, and some measurements were made 
through thin and broken clouds. An example is 
shown in Fig. 3. It shows an echo pulse from the 
ground, whose sharp edges aid in range determina- 
tion. It also shows time-resolved triple-echo pulses 
measured when viewing the ground over Ohio 
through two thin cloud layers. The first pulse in each 
triplet was the reflection from a cirrus cloud, while 


the second was the reflection from a midaltitude 
cloud and the third was from the ground. The differ- 
ent amounts of C0 2 column absorption for different 
path lengths are evident. Without a range-resolved 
receiver, the echo pulse signals and measurements 
from the three different path lengths would be 
mixed. Using the pulsed measurement approach 
allows using range gating in the data processing to 
isolate the surface-reflected signals from the preced- 
ing atmospheric backscatter, and hence eliminates 
path-length errors from scattering. 

After a series of test and checkout flights, four sci- 
ence flights were made in August. Diagrams of the 
flight’s ground-track patterns are shown in Fig. 4, 
and the flight locations, flight patterns, and altitude 
ranges are summarized in Table 2. Three of the four 
flights were coordinated with the NASA LaRC/ITT 
Acclaim C0 2 lidar team. They provide measure- 
ments from the AVOCET in situ sensor [29,30] and 
special radiosondes launched close to the time of 
the airborne C0 2 lidar overpasses. These flights 
permitted comparisons of measurements made at 
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Fig. 3. (Left) Typical laser pulse shape reflected from the ground as recorded by the photon-counting receiver. This was from accumu- 
lating 450 individual ground echoes integrated (histogrammed) over 0.9 s. (Right) Example of a raw (uncorrected) recorded backscatter 
profile when measuring to the ground through two thin cloud layers. It shows three time-displaced echo pulses per transmitted pulse. The 
smallest initial echo pulses are from the cirrus clouds just under the aircraft. The figure also shows the two range-displaced C0 2 column 
absorptions that are evident in the echoes from the midlevel cloud layer and the deepest absorption in the longest path to the ground. 
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Fig. 4. Maps of the ground-track patterns for the flights analyzed: 
(a) ARM-4, (b) Illinois, (c) North Carolina, and (d) Virginia. The 
flight locations and other parameters are summarized in Table 2. 


different altitudes and differential optical depths 
(DODs) to those calculated from in situ measure- 
ments and using the spectroscopic parameters for 
the C0 2 line from HITRAN 2008. 

5. Airborne C0 2 Measurements and Calculations 

In order to estimate the actual C0 2 column density 
during the flights, measurements of atmospheric 
temperature, moisture, and pressure vertical profiles 
were used from the radiosonde balloons, which were 
launched near the flight location. Their parameters 
were used in a 40-layer atmospheric model to com- 
pute dry air column density versus height to 
13 km altitude. The airborne flights were also coor- 
dinated with LaRC investigators, and the flights over 
the ARM site, North Carolina, and Virgina were 
underflown with in situ gas measurements made 
from the NASA LaRC UC-12 aircraft. This used 
the NASA LaRC atmospheric vertical observation 
of C0 2 in the Earth’s troposphere (AVOCET) gas 


analyzer to measure C0 2 concentrations. It sampled 
air and C0 2 concentrations every second from takeoff 
to 8 km, the UC-12’s upper flight altitude, and back 
to the ground. When the spiral was performed over 
an airport, the in situ measurements were obtained 
down to the ground via a “touch and go” of the airport 
runway. 

All AVOCET measurements were references to 
GPS latitude, longitude, and altitude measurements 
made on the UC-12. The spirals were timed to 
coincide with overflights of the GSFC C0 2 sounder. 
Vaisala RS-92 radiosondes were also launched dur- 
ing the Oklahoma, North Carolina, and Virgina 
flights to obtain the meteorological parameters 
needed for converting the C0 2 mixing ratios mea- 
sured by AVOCET to the path-integrated C0 2 col- 
umn density needed to compare with the lidar. The 
RS-92 radiosondes have pressure, temperature, 
and relative humidity absolute accuracies from 
1080 to 1100 hPa of 1 hPa, 0.5°C, and 5%, respec- 
tively. These radiosondes were launched from the 
ground at the center of the UC-12 spiral location, 
and when possible, they were launched at the com- 
pletion of the UC-12 spiral. The activities of the 
UC-12, Lear-25, and ground radiosonde personnel 
were coordinated to ensure the most coincident data 
acquisition for comparison of the in situ derived C0 2 
data and the lidar measurements. 

Based on these readings, the C0 2 line shapes, the 
C0 2 column number density, and the averaged two- 
way DOD for the 1572.335 nm line were computed 
versus flight altitude for each flight. These calcula- 
tions provided a reference comparison for the air- 
borne lidar measurements and showed how the 
C0 2 line shapes and depths should vary with flight 
altitude and range to the surface. A summary of the 
airborne concentration measurements, and the com- 
puted C0 2 column density, is plotted in Fig. 5. They 
show that the C0 2 mixing ratio was similar for three 
of the four flights. 

Over Illinois the AVOCET reading was not avail- 
able, so an analysis was used based on the nearest 
ground-based in situ sensor, which indicated a 
~30 ppm C0 2 reduction near the surface due the 
growing surrounding crop land. Figure 5 also shows 
the C0 2 number density plotted versus altitude, as 
a difference from the US standard atmosphere with 


Table 2. Summary of 2009 Science Flights 


Notation 

ARM-4 

Illinois 

North Carolina 

Virginia 

Location 

Lamont, OK 

Homer, IL 

Ahoski, NC 

Melfa, VA 

Flight dates (all 2009) 

August 4 

August 13 

August 17, am 

August 17, pm 

Flight pattern 

Bowtie 

Box 

Bowtie 

Racetrack 

Surface elevation (m) 

315 

178 

30 

15 

Surface pressure (mbar) 

974 

998 

1021 

1020 

Altitude change in flight 

Ascending 

Descending 

Ascending 

Ascending 

Primary altitude range (km) 

3.0-12.6 

12.6-5 

4.9-11.5 

4.9-11.5 

Altitude step size (km) 

1.5 

2.5 

1.5 

1.5 

Number of altitude steps 

10 (3 repeated) 

3 

5 

5 

Time span of data record (s) 

6300 

3400 

4500 

6000 
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Fig. 5. (Left) In situ C0 2 measurements from Avocet for the four flights in August 2009: ARM-4, Illinois, North Carolina, and Virginia. 
(Right) Computed C0 2 column density versus altitude above the surface for the four flights, minus that of the US standard atmosphere 
with a 385 ppm mixing ratio. The primary differences in the plots are caused by the different average column air densities caused by the 
different surface elevations of the flight locations. 


a constant 385 ppm mixing ratio. It clearly shows 
that the conditions of warmer dry air at lower pres- 
sure during the ARM and Illinois flights reduced 
their dry air density and hence their total column 
C0 2 compared to the later East Coast flights. The 
AVOCET instrument also sampled water vapor, 
and the column-integrated water vapor mixing ratio 
is plotted in Fig. 6 for all flights. In this plot the mea- 
surements for the Homer, Illinois flight were pro- 
vided by a nearby radiosonde. The results show 
that in the air columns measured in the flights, 
the water vapor mixing ratios were similar and var- 
ied smoothly from ~2% near the surface to ~0.5% for 
columns with upper endpoints near 12 km. 

For this analysis, we first performed a line fit, as 
described in the section below, based on all 20 meas- 
urement wavelengths. We then report the DOD on 
the fitted line where the on-line wavelength is set 
to the line peak, and the off-line wavelengths are 



Water Vapor Volume Mixing Ratio (%) 


Fig. 6. Plot of the path-integrated water vapor in the C0 2 meas- 
urement column for the four flights as a percentage of the dry air 
column. The water vapor was measured by the aircraft in situ 
sensor in the spiral segments of their flights, except for Homer, 
Illinois, where it was measured by a nearby radiosonde. 


in the wings. Specifically for these flights we 
computed and report 

DOD(p&) = OD(1 m ) 

- {OD(/l p £ - AX) + OD (A p k + AX)}/ 2, (1) 

where the OD values are measured on the fitted line 
shapes, and the off-line ODs are measured at offsets 
AX = 100 pm from the line peak, as shown in Fig. 7. 
For these choices of wavelengths, calculations from a 
layered atmospheric model shown in Fig. 8 show the 
DOD(pk) varies almost linearly with altitude. 

6. C0 2 Measurement Processing 

For the flights the lidar recorded the time- and 
wavelength-resolved laser backscatter with 0.9 s 
integration time. In subsequent analysis, the line 
shape measurements at each flight altitude step 
were averaged with 2-10 measurements/average. 
We used a C0 2 line retrieval approach, shown in 
Fig. 9, based on the Gauss-Newton method [31] to 
analyze each altitude averaged line shape. This ap- 
proach has sufficient free parameters to model and 
correct for instrument effects, to fit the resulting 
C0 2 line shapes, and to estimate the corresponding 
C0 2 column densities and optical depth at each alti- 
tude. The C0 2 retrieval algorithm yields an estimate 
of the mean C0 2 column density over the laser path 
length based on line absorption strength taken from 
the HITRAN 2008 [32] database. The input observa- 
tions were ratios of the photon counts in the surface 
echo signals at each wavelength after they were nor- 
malized by an estimate of transmitted pulse energy. 
The error covariance matrix for the observed signals 
was diagonal and equally weighted for all but the 
first three wavelengths. Examples of calculated and 
observed (retrieved) line shapes for flights over the 
ARM site and over North Carolina are shown 
in Fig. 7. 
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Fig. 7. Example plots from the line fitting and analysis approach. (Left) Examples of measured wavelength samples (dots) and fitted line 
shapes in transmission (solid curves) overlayed for averaged measurements above North Carolina from 3.1 to 11.2 km altitudes, which are 
color-coded. The wavelength axis is centered on the line peak. (Right) Computed line shapes in optical depth (OD) for the corresponding 
colors on left-hand plot. The three wavelength samples used for the 2009 flight analysis are indicated. 


The algorithm requires several other inputs. First 
are estimates of the vertical profile of temperature, 
pressure, and water vapor versus altitude. These 
are used to calculate the approximate C0 2 line 
shape, and were estimated from gridded meteorologi- 
cal fields from the Goddard Modeling and Assimila- 
tion Office for the locations at the time of the flights. 
They are used to calculate the wavelength-resolved 
absorption spectra of atmospheric C0 2 based upon 
HITRAN 2008 for each 1 km altitude bin. 

The range to the surface must be known to com- 
pute the volume density of atmospheric C0 2 from 
the measurement of DOD of the column. Our lidar 
processing does this by recording the trigger time 
of the laser pulse emissions and by computing an es- 
timate of the reflectance time of the echo pulse. We 
use a cross-correlation algorithm to optimize the 
range determination. Recent analysis [ 23 ] has shown 
this usually produces estimates with uncertainties of 
a few meters over most surfaces. 

The signal consisted of the detected photons in 
each pulse of the wavelength scan. The algorithm ad- 
justed the computed column average C0 2 density 
(and mixing ratio) to minimize the error between 
the computed line shape and that sampled by the li- 
dar. It does this also using an instrument model that 
allows degrees of freedom in several model variables. 
For these experiments, the lidar’s wavelength (i.e., 
wavelength per laser pulse number) was modeled 
as a quadratic function, and the three wavelength 
coefficients were solved for, using the ground wave- 
length calibration as a prior constraint. The lidar’s 
baseline response, or system transmission, was also 
allowed to vary with wavelength via a quadratic 
dependence with pulse number. This can change 
during flight due to the changing aircraft cabin 
temperature and pressure. Using this approach to 
simultaneously determine C0 2 density and instru- 
ment model parameters considerably improved the 
measurement accuracy by solving and correcting 


for slowly varying instrument offsets and nonlinear- 
ities. The result was that the residual difference be- 
tween the fitted C0 2 lines and lidar measurements 
on all flights was small, and the average rms differ- 
ence per wavelength sample varied between 0.7% 
and 1.7%. 

7. Expected Measurement Performance 

The error of the lidar measurement of optical depth 
and path-integrated number density depends on its 
signal and noise characteristics and the magnitude of 
bias errors. A simplified treatment for this approach 
[18], for a uniform atmospheric path and target 
at a fixed range R , illustrates some important 



Fig. 8. Computed DOD(pk) versus altitude for the 2009 flights, 
based on in situ readings of dry air and C0 2 concentrations and 
using HITRAN 2008 spectroscopy. The shapes of the DOD(pk) 
for all the locations are quite similar, and vary quadratically 
(nearly linearly) with range. A computation based on the US stan- 
dard atmosphere 385 ppm mixing ratio is also shown. Since it has 
a higher dry air density, it has a larger DOD at a given altitude. 
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Fig. 9. Algorithm used to process the 2009 airborne lidar measurements. 


dependencies. A complete analysis must account for 
many factors, including variability in the lidar 
parameters, atmospheric temperature and pressure, 
laser speckle, changing surface reflectivity and 
range, etc. [ 14 , 15 ]. 

This campaign was the first one in which this lidar 
was fully functional over the full range of measure- 
ment altitudes. During it, the lidar’s operational 
modes were still being optimized, and so factors such 
as laser power, detector gain settings, and receiver 
discriminator settings were varied during the flights, 
and from flight to flight. These operational settings 
changed the instrument responses, and hence 
changed the measurement offsets (or biases) for 
the flights. Hence in this paper we focus on analyzing 
the random errors in the measurements. 

For vertical measurements made with the on-line 
signal near the absorption peak, we can assess the 
primary impacts of random errors starting with 
the derivation given in Appendix A. The primary 
result is that the near surface number density and 
volume mixing ratio can be estimated from 


DODi (z s ,Z) 

fti(0) = X C o2n 2 (0) = ot\ p , (2) 

Ki 


where a 1 = l/{2[^ 1 (i on ,0) - 0 i(A o ff,Z-z s )]}, the range 
to the surface, R 1 = Z-z s , and the other terms are 
defined in Appendix A. Equation (2) is consistent 
with the results from the absorption calculations 
shown in Fig. 8 that were based on a layered atmos- 
pheric model. They show for A on = A pk and l of f well off 
peak that the values of DOD(pk) varied almost 
linearly with range R = Z - z s . 

The impact of errors in the airborne lidar measure- 
ments can be determined from Eq. (2). If the individ- 
ual measurement precisions of DOD and R are high, 
then the errors in the measurements, A, are small 


relative to their mean values. For this case, one 
can write the individual measurements as additive 
deviations (errors) around their mean values, so 
that 


rii ? 

-i'O- 




( ADOD \ 



DODj ? 

«dodJi+ '), 

V DODj ) 

and 


Rif 



(3) 


Since the deviations are small and have mean 
values equal to 0, the mean value for the C0 2 density 
can be approximated as 


rii = ai 


DODi 


(4) 


For these experiments the errors in the dry air 
densities are assumed to be much smaller than for 
C0 2 . Therefore relative error in the number density 
(2) is similarly 

(ffMffl -(f)' ,5) 

and the corresponding relative statistical error is 


g(A/i x ) _ AA APOD,) ^(ARQ 
ni V (DODO 2 (i?i) 2 ’ 

where rj 2 (x) is the variance of x and rr(x) is its stan- 
dard deviation. 
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With our airborne lidar, the random error in range 
was typically 2-5 m. The relative error in range 
varied from 0.1%, at 2 km altitude, to ~0.05% at 
10 km. These relative errors are much smaller than 
those for DOD, so the measurement errors in DOD 
dominate. For this case, the fraction error of the total 
column C0 2 number density can be written as 


from the finite SNR. At higher altitudes, the on-line 
signal varies as exp(-7?), and so the relative error at 
high altitudes increases as exp(i?/2). Since the rela- 
tive error has limits that grow at both lower and 
upper altitudes, there is an altitude band with mini- 
mum errors, which for these experiments was 
8-10 km. 


o(krii) <t(ADOD) 
n x DOD 


DOD ' (sNR off + SNR on ) ' 

(7) 


Since the airborne lidar scans the C0 2 line at a 
frequency f\ as and averages for an integration time 
T int , the total number of pulse measurements aver- 
aged is M = fi a s • T int . The average SNRs for the 
off- and on-line wavelengths are 


and 


SNR off = s/M 


N S (A off) 

/A/ /off) + 2N n 


SNR on = VM 


NjV on) 

Vnmtw; 


( 8 ) 

(9) 


where N n represents the detected background and 
dark counts accumulated over the laser pulsewidth 
for each laser firing, and N s represents the detected 
signal photons for that wavelength for each laser 
firing. 

The relative error varies significantly with 
range. As for most IPDA lidar, the C0 2 absorption 
for measurements at A on reduces the detected 
signal compared to l off , so that SNR on is usually 
the limiting error source. For this case Eq. (7) simpli- 
fies to 


<t(ADOD) 1 / VA^on) + 2A/ \ 

DOD DOD \ /MA/a on ) ) 


Since for the wavelengths chosen OD on » OD off , 


DOD 


OD on = 


^(0) • R 

ax 


(ID 


When the lidar measurement is also signal shot 
noise limited, Eq. ( 10 ) becomes 

ff(ADOD) = expCra/O) • R/2a{) 1 

DOD n/0) -R/ax y/MN s Q l off ) ’ 

Equation (12) summarizes the effect of varying 
range and off-line signal levels on the relative error. 
For these experiments, the maximum count rate of 
the photon-counting receiver was limited to prevent 
saturation. This limited the value of N s (A 0 ff) in 
Eq. (12) at lower altitudes, where the optical depth 
of the absorption is also small. Hence for low alti- 
tudes the 1/R term in Eq. (12) magnifies the errors 


8. Airborne Measurement Results and Discussion 

After data processing, the airborne measurements 
per site were summarized in a common format for 
all four flights, and were plotted versus flight 
location. The DOD(pk) values determined from 
the 0.9 s integration time were plotted versus the 
time of flight, along with the range estimates (R) 
to the surface. The same DOD(pk) values are plotted 
versus R to assess the linearity and scatter of the 
measurements versus altitude. Other analysis of 
the 2009 flights showed that the average signal 
levels follow predicted values [33], and that the 
altimetry measurements over water had an uncer- 
tainty of <3 m [23]. 

The 0.9 s measurements of DOD and range 
are plotted versus elapsed time for each flight in 
Figs. 10(a) , 12(a) , 14(a) , and 16(a) . Generally they 
show the range measurements and DOD(pk) varied 
together as expected on all flights, with the DOD(pk) 
and range being stable in the straight segments and 
both increasing due to the longer slant path ranges 
during the aircraft turns. The DOD(pk) readings 
have more scatter for measurements above 9 km. 
The DOD is plotted versus slant range for each of the 
flights in Figs. 10(b) , 12(b) , 14(b) , and 16(b) . The plots 
show the DOD generally followed the range as pre- 
dicted for all flights, and the DOD estimates have 
more scatter at ranges above 9 km. This is expected 
due to both the larger R 2 signal losses and the 
greater C0 2 on-line absorption for the longer 
columns. 

To determine statistics, for each flight the mea- 
surements were then segmented into groups of 
100 s time duration. For the recording period of 
2 s, these groups could have at most 50 individual 
(0.9 s duration) measurements. Screening criteria 
were then applied based on pulse energy, range, 
DOD, and others to remove outliers. Limiting the 
maximum spread in range worked well to remove 
measurements made in aircraft turns. Typically 
30-35 measurements survived screening in the 
groups. 

From these the mean values, standard deviations 
and relative errors were computed and plotted 
versus altitude. To allow assessing the equivalent 
measurement precision in an equivalent mixing 
ratio, the relative errors in the DOD values in 
Eq. (12) were scaled to equivalent errors in X C02 
by using 


*( AXc02) = *C02 • = Z C02 • ff( tS2 D) . (13) 


n X 


DOD 
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Fig. 10. Measurements from the flight above DOE ARM site, 
(a) Lidar measurements of range to surface (blue dots, left-hand 
axis) and DOD(pk) (red dots, right-hand axis) versus time for 
the same flight. The stair step changes in altitude are evident, 
and the DOD(pk) values follow the range until ~8 km, and then 
have more scatter above 8 km. A total of 1747 measurements (each 
1 s) were made, (b) Same lidar measurements of DOD(pk) plotted 
versus measured range to the surface. The line shown is a best 
quadratic fit to the measurements. The measurements show a 
good fit to the line, with more scatter above 10 km. 


This approach assumes the relative error in the 
measured C0 2 column absorption is much larger 
than that computed for the dry air column. Based 
on the airborne in situ measurements shown in Fig. 5, 
Xq 0 2 was taken to be 385 ppm. These relative errors 
were plotted versus range in Figs. 11(a) , 13(a) , 15(a) , 
and 17(a) . These plots show a significant dependence 
on the column length, as is expected from Eqs. ( 12 ) 
and (13). For ranges <8 km, the DOD is smaller 
and so the errors scale inversely with range. At 
ranges >10 km, the lidar has a declining signal 
strength range, due to both the R ~ 2 factor and on-line 
gas absorption, and so the random errors increase 
with altitude. As a consequence, in all flights the al- 
titudes with minimum relative errors were 8-10 km. 



of(AX C02 ) 



C7(AX C02 ) 

Fig. 11. Analysis of measurements from flight over ARM site, 
(a) Distribution of random errors calculated as 0(X CO2 ) versus 
altitude for measurements grouped every 100 s and for selected 
standard deviation ( R ) < 100 m. For low altitudes, the relative 
error increases because of a lower mean (DOD), while for upper 
altitudes the error increases due to declining received photons. 
The altitude with least scatter was 8-10 km, which had values 
of ~0.5 ppm. (b) Histogram of the same statistic for all altitudes. 
There were 38 groups, and the median error was 1.2 ppm. 


For this optimum altitude band most errors were be- 
tween 0.2 and 0.9 ppm. Histograms of relative errors 
for all ranges are plotted in Figs. 11(b) , 13(b) , 15(b) , 
and 17(b) . The lidar performance was similar on all 
flights, and some statistics from the analysis are 
summarized in Table 3. 

An analysis of the range measurements is also 
summarized in Table 3. It summarizes the range 
variability within the groups. This was primarily de- 
termined by the variability of the surface topography 
in the ~14 to 18 km distances traversed by the air- 
craft during the data accumulation time of 100 s 
for the groups. At higher altitudes the aircraft speed 
was higher and so the distances were longer. The 
flights over Illinois and North Carolina were similar 
with median range standard deviations of 16.7-19 m. 
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Fig. 12. Results from measurements from flight over Illinois, 
(a) Time history of lidar measurements of range to the surface (blue 
dots) and DOD(pk) (red dots). The increased range and increased 
DOD for the measurements during the turns (aircraft banks) at 
the corners of the box pattern are evident, (b) Plot of measured 
DOD(pk) versus measured range for the flight above Illinois. The 
line shown is a best quadratic fit to the measurements. 


Fig. 13. Analysis of measurements from flight over Illinois, (a) Al- 
titude dependence of (t(X C02 ) calculated for measurements. Mea- 
surements were evaluated in groups of 100 s elapsed time and 
screened for only those groups whose standard deviation (range) 
is <100 m. Typically each group has 32 measurements, (b) Histo- 
gram of the same statistic for all altitudes. There were 16 groups 
with a median error of 1.3 ppm and 11 groups with errors <1 ppm. 


The ARM-4 flight and the Virginia flight were also 
similar with median values of ^10 m. The ARM 
site had a minimum range standard deviation of 
3 m, while the Virginia flight, which was partially 
over the flat Chesapeake Bay water surface, had a 
minimum range standard deviation of 1 m. For the 
8-10 km altitudes where performance was best, 
the relative range variability was small, and the 
dominant error was measuring DOD(pk). 

9. Summary 

During the summer of 2009 we demonstrated air- 
borne lidar measurements of C0 2 absorption and 
column abundance using a pulsed direct-detection 
lidar based on the IPDA technique. Measurements 
were made at altitudes from 3 to 13 km over four 
locations and some through cirrus clouds. They 
showed clear absorption line shapes, which increased 
in optical depth with increasing aircraft altitude. The 


instrument's line shapes were estimated via a C0 2 
line shape retrieval algorithm, which permitted solv- 
ing and correcting for known instrument factors. The 
postprocessed line shapes agreed well with ones cal- 
culated from in situ measurements and radiosondes. 
The residual difference between the fitted C0 2 lines 
and lidar measurements was small, and the average 
rms difference per wavelength sample varied be- 
tween 0.7% and 1.7%. Analysis of the results on 
all flights showed good agreement with the nearly 
linear dependence of the DOD(pk) with altitude as 
predicted from calculations based on atmospheric 
models. The scatter in the relative error of DOD var- 
ied with aircraft altitude as expected, and the alti- 
tude range with lowest scatter was 8-10 km. Over 
the entire altitude range the median measurement 
precision varied from 0.9 to 1.2 ppm. In the optimum 
altitude range of 8-10 km, most measurements had 
precisions between 0.2 and 0.9 ppm. 
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Fig. 14. Measurements from NC flight, (a) Time history of the 
measured range to the surface (blue dots) and DOD(pk) (red dots). 
There were 1151 one-second measurements for this flight, (b) Plot 
of the measured DOD(pk) versus the measured range for the same 
flight. The line is a best quadratic fit to the measurements. 

After these flights we made additional flights on 
the NASA DC-8 during July 2010 to 13 km altitudes 
with some instrument modifications, which im- 
proved the measurement performance. For these 
campaigns, the lidar settings were better controlled, 
which allowed assessments of instrument bias levels. 
Those measurements have been analyzed [34] and 
will be reported in the future. 

Appendix A: Measurement Model for Random Errors 

Here we derive a simplified measurement model to 
allow evaluating the characteristics of random errors 
for these airborne measurements. The following 
discussion is for measurements using only two wave- 
lengths, while our line-scanning approach steps 20 
wavelengths across the line shape and uses a least- 
squares fit of a predicted line shape to the measure- 
ments. There are several advantages of the using 
additional measurement wavelengths distributed 
across the line. The additional information from the 


Fig. 15. Analysis for NC flight, (a) Altitude dependence of calcu- 
lated g(X c0 2 ) for measurements. Measurements were evaluated in 
groups of 100 s elapsed time and screened only for those groups 
whose standard deviation (range) is <100 m. Typically each group 
has ~32 measurements, (b) Histogram of the same statistic for 
measurements at all altitudes. There were 21 groups with a 
median value of 1.1 ppm, and 11 groups with relative error 
<1.1 ppm. 

line fit allows solving for spectroscopic shifts of the 
line with altitude, as well as instrumental offsets, 
wavelength dependencies, and nonlinearities. Hence 
the wavelength-distributed approach provides addi- 
tional information and a significant capability to solve 
for and thus limit bias errors. 

After applying corrections for the solved-for quan- 
tities, we have obtained similar results in analyzing 
random errors to those of the two-wavelength ap- 
proach if one segments its wavelength-scanned 
measurements into groups. Combining the signals 
from the four wavelengths closest to the line center 
as the on-line, and using the six the measurements 
nearest the wings (three on the short wavelength 
side, three on the long) gives similar results [35]. 

The volume mixing ratio of C0 2 in dry air can be 
written as 
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Fig. 16. Measurements from Virginia flight, (a) Time history of 
the measured range to the surface (blue dots) and DOD(pk) (red 
dots), (b) Plot of the individual 1 s measured DOD(pk), plotted 
versus the measured range for the same flight, along with a best 
quadratic fit to the measurements. 


Xq 02 ~ n l/ n 2, (Al) 

where n is the gas volume density, C0 2 is gas 1 and 
dry air is gas 2. For our case we have an IPDA lidar 
that measures the differential absorption in a round- 
trip pass of the laser beam through the column of 
C0 2 that extends from ranges 0 to R. For it the 
two-way transmission, t 2 (A, 0 ,R), at wavelength A on 
on the peak of the line is related to that at the wave- 
length in the wing of the line by 

T?(A 0 n, 0i? ) = 0,#) exp | -2 jf M-V.r) 

-^lOofr. r)] x n!(r)dr| , (A2) 

where ^i(A on ,r) and ^i(A 0 ff ,r) are the respective gas 
absorption cross sections at location r. 

For these experiments we have a nearly vertical 
nadir-directed path, so r = z. We also can set z = 0 
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Fig. 17. Analysis for the for the Virginia flight, (a) Altitude 
dependence of calculated g(X C02 ) for measurements. Measure- 
ments were evaluated in groups of 100 s elapsed time and screened 
only for those groups whose standard deviation (range) is <100 m. 
Typically each group has ~32 measurements. The relative errors 
are minimum at 8-10 km altitudes, (b) Histogram of the same 
statistic for the groups at all altitudes. There were 40 groups an- 
alyzed, the median error was 1.0 ppm, and 19 groups had <1 ppm 
error. 


to be the elevation of mean sea level (msl) and denote 
the elevation of the surface relative to msl as z s . For a 
vertical measurement from an aircraft at altitude Z 
to a surface at z s , the DOD of the C0 2 line can be 
written as 


DOD i(z s ,Z) = - In 


-£ 


( r\(Xo n ,z s ,Z) \ 

l^iOW ,2 s ,Z)j 


ki r) - <ti(A, off , r)] x m(z)dz 

(A3) 


= OD on (z s ,Z) - OD off (z s ,Z), 


(A4) 
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Table 3. Summary of 2009 Airborne Measurement Results 


Notation ARM-4 Illinois North Carolina Virginia 


Center of flight Location 

Number of 1 s measurements passing screening 
DOD(pk) versus range : linear fit results: R value of line fit 

Statistics for groups: 

Number of groups with std dev (range <100 m) 

Median number of measurements per group 
Median value of X Co2 statistical error (ppm) 

Altitude of min values (km) 

Minimum value of X Co2 statistical error (ppm) 

No. of groups with X Co2 statistical errors <1 ppm 
Std dev of range measurements in groups 
Median value (m) 

Minimum value (m) 


Lamont, OK 

Homer, IL 

Ahoski, NC 

Melfa, VA 

1747 

1182 

1151 

2031 

0.994 

0.989 

0.993 

0.993 

38 

16 

21 

40 

32 

32 

32 

32 

1.2 

0.9 

1.1 

1 

8-10 


8-10 

8-10 

0.6 

0.6 

-0.5 

-0.5 

17 

11 

11 

19 

10.1 

19.0 

16.7 

9.8 

3.0 

4.5 

5.1 

1.0 


where OD on and OD off are the optical depths of 
the gas absorption at the on and off-line points, 
respectively. 

Equation ( A4 ) was computed versus altitude using 
a layered atmospheric model for each flight location, 
and the results are shown in Fig. 8. The plot shows 
an almost linear dependency of DOD with altitude. 
To understand this dependency and assess errors, 
it is helpful to evaluate Eq. ( A3 ) by considering the 
terms separately. The gas absorption cross section 
can be written as 

o\ (A, r) = S 1 (r)/*! (A, r ) , (A5) 

where S(r ) is the absorption line strength. For the 
1572.33 nm line used, HITRAN shows that S(r ) in- 
creases for colder air, and so it increases modestly 
with higher altitudes. The line shape function f 1 
has unit area and depends on p(r) and T(r), the 
atmospheric pressure, and the temperature at 
location r. For measurements in the Earth’s atmos- 
phere, fi can usually be represented as a Voigt line 
shape. This is a convolution of the Lorentz line shape, 
caused by pressure broadening, which dominates 
in the lower atmosphere, and the Doppler line shape, 
caused by thermal motion of the gas mole- 
cules, which dominates in low pressures at higher 
altitudes. 

For these flights, Fig. 6 shows the C0 2 mixing ratio 
was nearly constant in vertical segments of the 
atmosphere sampled, so that 


n 1 ^X C02 n 2 . (A6) 


From the ideal gas law, the dry air density varies 
with altitude as 


n 2 (z) = 7i 2 (0) 


Pd(z) T (0) 
P d (0) T(z) 9 


(A7) 


where p d is the partial pressure of dry air. The total 
atmospheric pressure p = p d + p w , where p w is the 
partial pressure of water vapor. This was calculated 


from in situ measurements during the flights, and 
the results are plotted in Fig. 6. They show the 
column average of p w varied from 0.02 *p d for the 
lower to 0.005 * p d for the higher flight altitudes; 
hence for the altitude range of these flights p ~p d . 

As p(z) decreases with altitude, /q narrows. Since 
it has unit area, near the absorption peak, the line 
shape function value varies inversely with linewidth, 
so that /■(A 0 n) ex ( T(z)) n /p d (z ). Here n is the temper- 
ature exponent for the pressure-broadened line- 
width, which the HITRAN database gives as 
n = 0.67 for our C0 2 line. The result is that the 
change in line shape near the peak approximately 
compensates for n 2 (z) decreasing with altitude. So, 
for this special case, the on-line optical depth 


OD on (z s ,Z) * 2 


L 




T(zy 


T( 0)J 


which can be simplified to 


(A8) 


OD on (z s ,Z)K2X C o 2 n 2 (0) P S' 1 (z)f i (A pk ,0)dz, (A9) 

Jz s 

with a temperature-modified line strength 


S\( 2 ) = S 1 (z) 


T(z )■ 


T(0). 


72-1 


(A10) 


Using these results, Eq. ( A8 ) can be approximated 
as 


OD on (z s ,Z) ~ 2 (Z - z s ) c> i(A p k, 0)XcQ 2 n 2 (0). (All) 

Here (A pk , 0) = (z)fi (A pk , 0) denotes that the 

temperature-modified cross section is evaluated at 
the path-averaged temperature but at altitude equal 
to 0. Although Eq. ( All ) is useful for this error analy- 
sis, the linear relationship is approximate, due to the 
increasing value of S'^z) with z that causes a slightly 
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increasing slope of OD with z. This slight curvature is 
evident in the plots shown in Fig. 8. 

The optical depth of the off-line term is 

OD off (z s ,Z) = 2 f s 1 (z)/‘ i (; l off ,z) 

Jz s 

xX CO2 n 2 (0)^^dz. (A12) 

As shown in Fig. 7, for the A off chosen, 
OD off OD on . Computations show its value varies 
from 5% of OD on near the surface to <1% OD on at 
the upper flight altitudes. Given its minor role in 
Eq. (A4), we can replace its cross section by its 
path-averaged value: 


OD off (z s ,Z) « 2(Z - z s )cji(A 0 ff,Z - z s )X CO2 n 2 (0). 

(A13) 

Substituting Eqs. ( All ) and ( A14 ) into Eq. ( A4 ) 
yields 

DOD^Z) = 2(Z - z s )[di(kp k , 0) 

- ^(i off ,Z-z s )]X CO 2^(0). (A14) 

Hence for this case, for a column with approxi- 
mately a uniform mixing ratio, the DOD also changes 
approximately linearly with the range to the surface. 
This result is consistent with the computed values 
shown in Fig. 8. 

From Eq. (A14), the C0 2 volume density at the 
surface can be determined from measurements of 
DOD and range by 

DODxfe^Z) . 

^i(O) = ^co 2 ^ 2 ( 0 ) = a 1 — , (A15) 

K 1 


where a 1 = l/{2[er 1 (A on , 0) - tfi(A 0 ff,Z - 2 S )]} and the 
range to the surface R 1 = Z -z s . 
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